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From the demonstration of saturable absorption by Vavilow and Levshin in 1926 [l|, and with 
invention of the laser, unavailability of strongly nonlinear materials was a key obstacle for develop¬ 
ing optical signal processing, in particular in transparent telecommunication networks. Today, most 
advanced photonic switching materials exploit gain dynamics [2] and near-band and excitonic effects 
in semiconductors [3], nonlinearities in organic media with weakly-localized electrons [4] and nonlin¬ 
earities enhanced by hybridization with metamaterials [5]. Here we report on a new type of artificial 
nonlinearity that is nano-opto-mechanical in nature. It was observed in an artificial metamaterial 
array of plasmonic meta-molecules supported by a flexible nano-membrane. Here nonlinearity is 
underpinned by the reversible reconfiguration of its structure induced by light. In a him of only 100 
nanometres thickness we demonstrated modulation of light with light using milliwatt power telecom 
diode lasers. 



FIG. 1: Nonlinearity in nano-opto-mechanical meta¬ 
material. Light-induced (red) forces in the plasmonic meta¬ 
material array cause nanoscale reversible displacements of its 
small and light building blocks supported by an elastic nano¬ 
membrane. They move fast, with a response time reaching mi¬ 
croseconds. These displacements change the plasmonic spec¬ 
tra of the metamaterial array and its transmission. This is 
used to modulate another, weaker beam of light (green) at a 
different wavelength. 

Some exceptional opportunities for developing en¬ 
gineered nonlinear media are provided by nano- 
optomechanics that take advantage of the changing bal¬ 
ance of forces at the nanoscale |5ffi3]. With the decrease 
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in the physical dimensions of a system the electromag¬ 
netic forces between constituent elements grow, as may 
be illustrated by the repulsion of electrons as their sepa¬ 
ration diminishes. In contrast, elastic forces, such as the 
force restoring a compressed spring, decrease with size. 
Moreover, the nanoscale metamaterial building blocks 
have high natural frequencies and thus can be moved 
very fast, potentially offering GHz switching bandwidth 
for elements of sub-micron size. 

Recently this was exploited to develop reconfigurable 
metamaterials fabricated on nanoscale elastic mem¬ 
branes. These photonic metamaterials can be driven 
thermo-elastically [14j and with electromagnetic forces, 
such as the Coulomb force between charged elements of 
the nanostructure m, or the Lorentz force acting on 
currents running through conductive elements in mag¬ 
netic field [1.6]. Their electro-optical and magneto-optical 
switching characteristics surpass those of natural media 
by orders of magnitude. 

Metamaterials are in essence arrays of optical res¬ 
onators where light-induced nano-opto-mechanical phe¬ 
nomena can play crucial roles. For instance it has been 
suggested that nearfield PH and Casimir forces mm 
in metamaterials are modified in the presence of light. 
Forces between oscillating plasmonic or displacement 
currents induced by light in metamolecules were also the¬ 
oretically shown to be sufficient to drive reconfiguration 
of a metamaterial structure in the optical part of the 
spectrum [2IH22], and proof-of-principle experiments on 
such interactions between individual metamolecules have 
been reported at microwave frequencies [23]. 

Here we show that reversible reconfiguration of a 
plasmonic metamaterial nanostructure driven by opti¬ 
cal forces between its illuminated elements can be the 
source of a strong optical nonlinearity. We show that 
this cubic nonlinearity may be used to modulate light 
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with light at milliwatt power levels. Moreover, although 
for the majority of media the magnitude of nonlinearity 
tends to be proportional to its response time, the nano- 
opto-mechnical nonlinearity is three orders of magnitude 
faster than could be expected from this otherwise univer¬ 
sal trend. 

The plasmonic nanomechanical metamaterial was fab¬ 
ricated on strips of suspended dielectric membrane of 
nanoscale thickness in such a way that plasmonic el¬ 
ements of the metamolecules were located on differ¬ 
ent strips. Here light illumination leads to electromag¬ 
netic and thermal forces that induce nanoscale strip dis¬ 
placements and thus reconfiguration of individual meta¬ 
molecules in a way that affects their plasmonic reso¬ 
nances, see Fig. [l] This leads to modulation of the 
metamaterial’s transmissivity at different wavelengths. 
The nonlinear response of metamaterial, that is only 100 
nanometres thick, can be observed at only a few milli¬ 
watts of continuous laser power, using lasers operating 
at telecommunication wavelengths. 

We developed nano-opto-mechanical metamaterial on 
the basis of a n-shaped resonator design known for ex¬ 
hibiting plasmon-induced transparency [24l|25], see Fig.Jl] 
and Fig. UK To allow mechanical deformation of the 
700 nm x 700 nm plasmonic n metamolecules, the hori¬ 
zontal and vertical gold bars were supported by different 
flexible silicon nitride strips of 28 fi m length spaced by 
alternating gaps of 95 nm and 145 nm. The nanostruc¬ 
ture was fabricated by focused ion beam milling from a 
50 nm thick silicon nitride membrane covered by 50 nm 
of gold. Such structures can be reconfigured by (i) dif¬ 
ferential thermal expansion through light-induce heating 
of the bimorph gold and silicon nitride layers and (ii) by 
near-field optical forces acting between elements of illu¬ 
minated plasmonic resonators. 

Experimental measurements and full 3D Maxwell sim¬ 
ulations show a pronounced near-infrared absorption res¬ 
onance around 1240 nm, Fig. and Fig. |Jd. Maxwell 
stress tensor calculations [20] assuming normal incidence 
illumination of the non-diffracting, periodic metamaterial 
reveal optical forces acting on the n-resonators around 
these resonances, see Fig. Upc. As the normally incident 
photons only carry momentum along the z-directionthe 
net force along z must comply with the momentum 
transfer associated with absorption A and reflection i?, 
F z i + F z 2 = (A + 2R)P/c , where P is the incident power 
per unit cell area and c is the speed of light in vacuum. 
In close agreement with these relationships, our simula¬ 
tions show substantial differential optical forces F 2 — F\ 
between the unit cell’s strip segments reaching about 
0.4 P/c along y and 2.8 P/c along z. 

Optical forces in such metamaterial nanostructures can 
be understood as time-averaged Coulomb and Lorentz 
forces acting between the currents and oscillating dipole 
charges of the plasmonic resonators in the presence of 
the illuminating electromagnetic wave [22]. For example, 
repulsive and attractive optical forces within the meta¬ 
material plane for illumination at wavelengths of 1310 nm 




FIG. 2: Optically reconfigurable plasmonic meta¬ 

material. (a) Scanning electron microscope image of 
the metamaterial nanostructure consisting of asymmetrically 
spaced gold (yellow) resonators on silicon nitride strips (red). 
(b) Transmission T, reflection R and absorption A spectra of 
the metamaterial for x-polarized light incident on the silicon 
nitride side. Dashed lines indicate the optical 1310 nm probe 
and 1550 nm pump wavelengths used in the modulation mea¬ 
surements of the main manuscript. The inset shows detailed 
unit cell dimensions in nanometers. The scanning electron 
micrograph was recorded by FEI Helios 600 NanoLab and 
the optical spectra were measured using a microspectropho¬ 
tometer (CRAIC Technologies). 


and 1550 nm correspond to repulsive and attractive in¬ 
teraction of the dipole charges on neighboring strips, see 
Fig- [3}i ,b. Forces normal to the metamaterial plane result 
from a combination of optical pressure and the Lorentz 
force acting on the moving charges of the plasmonic mode 
in the presence of the incident wave’s magnetic field. 

We investigated the nonlinear optical properties of the 
metamaterial in a pump-probe experiment as presented 
in Fig. UK As the pump and probe optical sources we 
used fiber-coupled telecommunication laser diodes oper¬ 
ating at the wavelengths of 1550 nm and 1310 nm, re¬ 
spectively. The intensity of the pump beam was modu¬ 
lated by a fibre-coupled electro-optical modulator. The 
pump and probe were then combined into a single beam 
using a fiber coupler, de-coupled into free space and fo¬ 
cused on the sample placed in a microscope using fo¬ 
cusing and collection objectives. The sample was placed 
in a cell evacuated to the pressure of 30 Pa to prevent 
damping of the nano-mechanical motion. Light trans¬ 
mitted through the sample was detected by an InGaAs 
photodetector and a lock-in amplifier locked to the pump 
modulation frequency. In order to optimize the light- 
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FIG. 3: Optical spectra and optical forces, (a) Simu¬ 
lated metamaterial transmission T, reflection R and absorp¬ 
tion A spectra. Insets show maps of the instantaneous values 
of the optically induced charges at the probe and pump wave¬ 
lengths. (b) In-plane-of-metamaterial component of optical 
forces F y acting between the strip segments of an individ¬ 
ual unit cell according to Maxwell stress tensor calculations, 
(c) Normal to the metamaterial component of optical forces 
F z 1,2 acting on the strip segments along the light propagation 
direction. The total optical force on the unit cell (green line) 
is presented alongside the value expected from reflection and 
absorption (dotted black line). Forces are shown per unit cell 
in units of P/c, where P is the incident power per unit cell 
and c is the speed of light in vacuum. All quantities are shown 
for x-polarized illumination of the silicon nitride side. 


induced nano-opto-mechanical modulation, we pumped 
the nanostructure at 1550 nm just above its absorption 
resonance, where simulations predict significant in-plane 
forces. 

At low modulation frequencies, below 100 kHz, the 
optical pump leads to pronounced modulation of the 
metamaterial’s transmission at the probe wavelength, 
see Fig. [4^. For 0.9 mW pump power (peak intensity 
920 W/cm 2 ) a modulation amplitude on the order of 1% 
is observed at 25 kHz modulation. The low-frequency 
component of the nonlinear response has a thermal na¬ 
ture and drops rapidly with increasing modulation fre¬ 
quency, fading at a few 100s of kHz. In this frequency 
range, the modulation amplitude halves as the modula¬ 
tion frequency doubles which is consistent with a ther¬ 
mal deformation of the bimorph Au/SiNx membrane, see 
inset to Fig. ^p. This is consistent with thermal cal¬ 
culations for 50 nm gold f26i and silicon nitride m 
layers, which predict conductive cooling timescales on 
the order of 20-30 /is. The light-induced heating of 
the nanostructure is balanced by heat conduction along 
the metamaterial strips and according to our calcula¬ 
tions can reach hundreds of degrees with pump excita¬ 
tion at 1 kW/cm 2 . As the thermal expansion coefficient 
of gold (14.4 xlO 6 /K) exceeds that of silicon nitride 
(2.8 x10 _6 /K) 5-fold such large temperature changes 
will bend the strips. When heated by the modulated 
pump laser narrow strips will rise higher than wider strips 
due to more complete gold coverage; the metamolecule 
changes its shape and therefore its plamonic absorption 
and transmission for the probe wavelength change. 

While the thermal nonlinearity easily dominates at low 
modulation frequencies, other modulation mechanisms 
are at play at higher frequencies. Indeed, as differential 
thermal expansion in the layered structure can only drive 
motion out of the metamaterial plane, it cannot drive 
in-plane resonant oscillations of the metamaterial strips. 
Several resonant peaks of nonlinear response are seen at 
much higher frequencies around 600 kHz, 1 MHz and 
2 MHz, see Fig. [ 4 J 3 . Mechanical eigenmode calculations 
identify the resonances as the fundamental modes of os¬ 
cillation of the narrow and wide strips, corresponding to 
vibration normal to the metamaterial plane (around 600 
kHz) and motion within the metamaterial plane (at MHz 
frequencies), see insets. Comparing the sub-/is heat¬ 
ing and cooling cycles to the characteristic strip cooling 
timescale, we argue for their non-thermal origin. 

To assess the origin of the high frequency modulation 
we estimate the optomechanical deformation that can be 
expected from the optical forces acting between different 
parts of the metamolacular resonator. Based on the re¬ 
sults of simulations presented in Fig. Hf pump-induced 
optical forces reach 0.5 P/c (0.2 P/c ) for the narrow 
(wide) strips. For a strip pair with 29 unit cells and 
920 W/cm 2 pump intensity, this corresponds to about 
220 fN (90 fN) force per strip. Estimating the spring 
constant k from strip mass m and measured resonance 
frequency / as k = (2 tt f) 2 m =28 fN/pm (56 fN/pm) 
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FIG. 4: Modulating light with light using nano-opto-mechanical metamaterial, (a) Schematic of the pump-probe 
experimental setup, (b) Transmission modulation depth at A = 1310 nm as a function of pump power for several pump 
modulation frequencies and pump modulation depth close to 100%. The inset shows a hyperbolic fit ~ 1 // (black dashed line) 
of the thermal modulation tail for a pump power of 0.9 mW. (c) Modulation depth as a function of pump laser power and 
modulation frequency. Corresponding mechanical modes of the metamaterial nanostructure are shown near resonant peaks. 


for out-of-plane deformation, we may expect the gap be¬ 
tween strips and thus between elements of the meta¬ 
molecule to change by 6 pm. Corresponding estimates 
for in-plane deformation predict a gap change on the or¬ 
der of 1 pm. Such modulation is not sufficient to create a 
detectable level of non-resonant high-frequency modula¬ 
tion. However, much larger displacements are expected 
at the structure’s mechanical resonances, where the dis¬ 
placement will be enhanced by the quality factor of the 
mechanical resonator. Our optical measurements show 
mechanical resonance quality factors on the order of 100, 
which may serve as a lower limit for the quality factor 
of mechanical resonances of individual strips as the over¬ 
all resonance profile is affected by inhomogeneous broad¬ 
ening due to slight structural variations from strip to 
strip. Taking this enhancement into account, optical 
forces are definitely sufficient to induce resonant strip 
displacements on the order of 1 nm, which is consistent 
with the observed transmission modulation on the order 


of 1%. Therefore we argue that the resonant response at 
600 kHz, 1 MHz and 2 MHz is predominantly driven by 
non-thermal, optical forces. 

Absorption in a nonlinear medium is conventionally 
described by the expression — ^ = al + /3I 2 + ..., where 
I is light intensity, z is the propagation distance and 
a and (3 are the linear and nonlinear absorption coeffi¬ 
cients. As the observed nonlinear transmission change 
is proportional to the pump power, we can quantify the 
nonlinearity of the reconfigurable photonic metamaterial 
by estimating its first nonlinear absorption coefficient (3. 
Assuming that the nonlinear transmission change AT re¬ 
sults from nonlinear absorption, f3 ~ A T/(It), where t is 
the metamaterial’s thickness. At 1 MHz modulation fre¬ 
quency, (3 ~ 5 • 10“ 3 m/W corresponding to a nonlinear 
susceptibility on the order of lm(x^)/n 2 ~ 10 -12 m 2 V -2 
and a nonlinear imaginary refractive index change Ak ~ 
5 • 10 -3 at 920 W/cm 2 pump intensity. We argue here 
that the observed nano-opto-mechanical nonlinearity is 
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FIG. 5: Optical nonlinearities in metamaterials and 
conventional media. Cubic optical nonlinearities due to 
thermal effects, saturated absorption, molecular orientation 
and electronic anharmonicity follow a general trend of de¬ 
creasing achievable modulation frequency v with increasing 
magnitude of the nonlinearity: x!^ v ~ 10 -9 m 1 2 V -2 s -1 [28] . 
Metamaterial patterning of plasmonic metal (points 1 and 2, 
29] 30 ), hybridization of carbon nanotubes (point 3, [STl), 
graphene (point 4, [32]), and semiconductors (points 5 and 6, 
OH EU) with metamaterials resonantly enhances the magni¬ 
tude of the nonlinear response without affecting its response 
time. Nano-opto-mechanical metamaterial demonstrated in 
the present work also departs from the common trend ex¬ 
hibiting that is three orders of magnitude bigger than 

in conventional nonlinear media (point 7). 


special as it exhibits a departure from the prevailing rela¬ 
tion between speed and magnitude of nonlinear responses 
(Fig# This happens in a way that allows larger non¬ 
linearity or higher modulation frequency to be achieved 
than could be expected from traditional media. Indeed, 
it appears that for many mechanisms of optical nonlin¬ 
earity the achievable modulation frequency v is related 
to the magnitude of the nonlinerity x^ in such a way 
that x!^ v ~ 10 -9 m 2 V -2 s -1 [28] . At the resonance of 
in-plane motion, MHz rate modulation can be achieved 
at about 1 kW/cm 2 of light intensity corresponding to 
about x!^ v ~ 10 -6 m 2 V -2 s -1 . 

In summary, we provide the first experimental demon¬ 
stration of a giant nano-opto-mechanical nonlinearity in 
plasmonic metamaterial that allows modulation of light 
with light at MHz frequencies and miliwatt power lev¬ 
els. Bimorph deformation resulting from optical heating 
at low modulation frequencies and electromagnetic forces 
between elements of plasmonic resonators drive the non¬ 
linear response. The magnitude and spectral position of 
the underlying optical resonance can be controlled across 
the optical spectral range by design of the plasmonic res¬ 
onators and the supporting nano-membrane strips. 


I. METHODS 


A. Plasmonmechanical metamaterial fabrication. 

Starting with a commercially available 50 nm thick 
low stress silicon nitride membrane, a 50 nm thick gold 
layer for the plasmonic metamaterial was thermally evap¬ 
orated. The gold-coated membrane was structured with 
a focussed ion beam system (FEI Helios 600 NanoLab). 
The pattern of plasmonic resonators was milled and then 
the supporting silicon nitride membrane was cut into sus¬ 
pended strips with tapered ends. Detailed dimensions of 
the nanostructure are given by Fig. [2 Jd. 


B. Experimental characterization. 


The cubic nano-opto-mechanical nonlinearity was mea¬ 
sured by detecting the pump-induced modulation of the 
metamaterial’s transmission (Fig. [4^). The probe beam 
at wavelength 1310 nm was generated by a fiber-coupled 
telecom laser diode (Thorlabs FPL1053S). The pump 
beam from a fiber-coupled laser diode at wavelength 
1550 nm (Thorlabs FPL1009S) was modulated with a 
fiberized electro-optical modulator (EOSpace AX-0K5- 
10-PFA-UL) and then combined with the probe beam in 
a fiber coupler (JDSU L2SWM1310/1550BB). The trans¬ 
mitted beam was filtered to remove the pump and the 
probe signal was detected by an InGaAs photodetector 
(New Focus 1811) and a lock-in amplifier (Stanford Re¬ 
search SR844). 

In all optical experiments and simulations the electric 
field was x-polarized (parallel to the strips) and incident 
on the silicon nitride side of the sample. 
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